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Abstract 

It is shown that difficulties in reconciling the values of the cross sec- 
tion cr(e^e~ J /ip + rjc) measured at BELLE and calculated within non- 
relativistic QCD (NRQCD) are caused not by some misinterpretation of the 
data or other exotic explanations, but by poor applicability of NRQCD to 
such processes. 

We use general theory of hard exclusive processes in QCD together with 
more realistic models of charmonium wave functions, and show that the 
BELLE result can be naturally explained. 
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1. A surprisingly large rate for hard exclusive processes of the type 
e'^e~ — > J/ip + T]c observed at BELLE P still remains unexplained. In 
this experiment the process e'^e~ —>■ J/ip + X was studied. The cross-section 
of e~^e~ —>■ J/ijj + rjc was then extracted from the number of events in the r]c 
peak in the mass spectrum of the system recoiling against the reconstructed 

In a recent upgrade of the BELLE analysis with a data sample of 155 fb~^ 
(21, [3 the cross-section of the e+e^ J/ip + i]c process has been found 
equal to (25.6 ± 2.8 ± 3A){h/Br{ric > 2 charged). BELLE also performed 
simultaneous fits to the production and helicity angle distributions. The 
measured angular and helicity distributions for J/ip + rjc have the general 
form (1 + acos^6') and are consistent with the expectations for production 
of this final state via a single virtual photon, aprod = Oihei = +1- 

From the theoretical side, this cross section was calculated in j3] within 
the NRQCD approach, and much smaller value ~ 2.3 fb was obtained. This 
large discrepancy initiated further studies, both experimental and theoretical. 
Various explanations were proposed, e.g. that the two-photon production of 
{J/ip + J/iJj) can be significant and can imitate those of J/ ip + rjc jH], or 
even more exotic variants of the scalar gluonium production (for its mass 
happily coinciding with the charmonium energy levels) etc. 

The main purpose of this paper is to show that the origin of the discrep- 
ancy is due to a poor accuracy of NRQCD when applied to such processes. 
The main reason is that the charm quark is not sufficiently heavy and, as 
a result, the charmonium wave functions are not sufficiently narrow for a 
reasonable application of NRQCD to the description of charmonium produc- 
tion. And, as usual, hard exclusive processes are particularly sensitive to the 
widths of hadron wave functions. Below we describe the properties of the 
model wave functions of charmonium (which we consider as more realistic 
in comparison with the extreme 5-function-like wave functions of NRQCD), 
and show that the value of a{e^e^ —>■ J/tp + rjc) measured at BELLE can be 
naturally obtained. 

2. The cross section of the process e+e^ 7* ^ J/'ip^iPi) + Vc{P2) has 
the standard form: 




^ Later, this possibility was excluded by BELLE j2|. 
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where {\p\/EY is the P-wave phase space factor and Qc = 2/3 is the charm 
quark charge, while the form factor is defined as: 

{ip±{pi),vM\JM\^) = e^.p.e'ip1p^,Fvp{s) . (2) 

Since there is one form factor only, the angular distribution is pure kinemat- 
ical: ~ (1 + cos^ 6). 

We will need only the asymptotic form of Fvp{s) because s = M^^^^^^ ~ 
112 GeV^ in the BELLE experiment. 

General theory of hard exclusive processes in QCD has been developed 
in |Zj;|H] (see P for a review). In |7j it was obtained that at large s = 
{Pi +^2)^ the leading power term of the general two-hadron form factor has 
the following behaviour (up to logarithmic corrections): 

(i7i(pi,si,Ai; H2ip2,S2,X2)\Jx\0) ~ i^-j=j . (3) 

Here: Hi and H2 are any two hadrons with momenta pi^2 , spins Si^2 and 
helicities Ai^2 in c.m.s., the current helicity (here, the projection of the photon 
spin onto direction pi) is: A = (Ai — A2) = 0, ±1. n.mm is the minimal number 
of point-like constituents (quarks or gluons) in these hadrons, n^m = 2 for 
mesons and n^m = 3 for baryons. It is seen that the asymptotic behaviour 
is independent of hadron spins, but depends essentially on their helicities. 

For the process considered, e~^e~ ^ 7* — J/4'±{pi) +^7c(P2), the J/ip- 
meson is transversely polarized, i.e. has helicities |Ai| = 1 only. So, the 
matrix element in eq.(3) behaves as ~ 1/s. Since in eq. (2) e± ~ 1, while 
Pi ~ P2 ~ all that results in Fvp{s) ~ 1/s^. 

The leading term of Fvp{s) is given by four similar diagrams, one of which 
is shown in Fig.l. Its explicit form will be given below in eq.(9,10), but we 
describe first the properties of the meson wave functions entering eq.(lO). 

3. The twist 2 and twist 3 light cone wave functions Vi{x) and Pi{y) of the 
^Si and ^Sq states of quarkonium made of the heavy QQ quarks are defined 
in the standard way (see e.g. [H]), we neglect higher twist wave functions 
giving only power corrections in comparison with those written below. For 
simplicity, we also do not distinguish My from Mp and use M which can be 
taken as an appropriate average, for instance: M = {SMy + Mp)/4. 
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Figure 1: One of the four similar diagrams for the form factor Fjj^j^^. 
For the meson with the hehcity A: ^ 

+ P^^V^:,) + ft^^*)^^ e^^p'^)VT{x) + a^^^^^^ . (4) 

Taking A = in eq. (4) and using Me^^o — > for large p, one obtains the 
standard definition of the leading twist wave function Vl{x) of longitudinally 
polarized vector meson jH]: 

{Vx=o{p)\Ciz)^,C{-z)\0) = fv C dx, e^(P^)(^^-^^)\4(x), (5) 

J o 

so that: V{x) = - V^{x). 

Xi and X2 = (1 — Xi) in eqs.(4,5) are fractions of the meson momentum 
p~^ = {E + Pz) = qo carried by quarks at large p^, i.e. [+ , ±, — ] components 
of the quark momentum are: ki = [xiqo, 0±, Mq/xiq^ (and similarly for 
three other quarks), and we neglect the quark transverse momentum inside 
the heavy quarkonium in comparison with its mass. 

^ We use the conventions: 

O-^x,. = [7ai,7i^]/2, 75 = {l5)Bj-Dr, P = Pf,!^ = Polo- PI- 
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For the P-meson: 

J a/3 

The wave functions and operators entering eqs.(4-6) are defined at the soft 
scale /i*, /i* <^ qb, where qb = Cpasifi — qB)MQ/2 is the Bohr momentum 
of the heavy quark in the quarkonium. 

All wave functions in eq.(4-6) are symmetric under: xi ^ X2, yi y2 
and normalized: Jq dxiVi{x) = I, i =J-, L, T, A; J^dyiPi{y) = 1, i = 
A,P,T. 

The values of dimensionless constants in the above formulae follow di- 
rectly from the exact QCD equations of motion: iDQ = MqQ (see e.g. ch.9 
and Appendix C in [5]; here Mq is the "soft" mass of the heavy quark, e.g. 
the appropriately defined perturbative pole mass) ^ : 

/>-) = 5j^. /;(M-) = (Tf -A.). (7) 

However, the appropriate scale n for the wave functions entering the 
hard form factor Fyp{s) will be /i^ ~ k'^ rather than /x*^, see Fig.l. For 
simplicity, we neglect in what follows the complicated but slow logarithmic 
evolution of (normalized) wave function forms, and will account only for the 
overall renormalization factors of the local tensor and pseudoscalar currents 
and for running of the quark mass. Unlike the larger soft mass Mq, the mass 
Mq entering below the perturbative logarithmic renormalization factors and 
couplings fl^pik'^) will be the "hard" mass, e.g. Mq = M^'^{fi = M^^), 

^ As in eq.(7) originates from the matrix element: {V\{p)\Q{iD^ — «-D^)Q|0)^. = 
AsfvM^e^. It is 0{v^) parametrically, and up to 0(i;4): ~ (4M*2 - M^)/6MM*. 
For the c-quark pole mass M* = 1.65 GeV : Ag ~ 0.05, and will be neglected in what 
follows. 
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which is smaller because a part of the heavy quark self-energy in the interval 
(1/Mq) < R < (l//i*) is now excluded from Mq. So, we use in eq.(lO) (the 
renormalization factors Zp, Zt and Z^ are given below): ^ 



4. The leading contribution to the form factor Fvp{s) is calculated in 
a standard way (see ch.9 in jHj where the calculation of the form factor 
7* — i> + vr is described in detail ^ ) , and the result for the leading term 
has the form: 



32n fvfpM 

\Fvp{s)\ = — lo, (9) 

'io 



JO [ a{x,yjs[xj M d[x,yj s{xj 



l VL{x)PA{y) 1 {1 - 2y^) V^jx) PAjy) 

2 d(x,2/) 2 s{y) d{x,y) ^ ' 



Here: s = AE^, ql = {\p\ + Ef ~ (s - 2M\ = (A;i + h) is the gluon 
momentum in Fig.l, d{x,y) and s(a;) originate from the gluon and quark 
propagators, Zt and Zp are the renormalization factors of the local tensor 
and pseudoscalar currents: 

^(^,^) = |=(^i + ^)(2/i + ^), ^=(^^^1, (11) 




^ The function Pt{x) gives no contribution to the form factor Fvp{s). 

^ The calculation for charmonium is even simpler because the three-particle ccg- 
wave functions of twist 3 containing the additional valence gluon have couplings /J ~ 
jvl^QCol^i i-^- ^'^^ suppressed by a factor ^ {Kqcd/M)'^ in comparison to /yM ~ fpM 
in eqs.(4-6), and can be neglected. 



6 




Mgifx') = Mq, Zt = = A;2), Zl^ = Z^l^^ = a') , (13) 

where Mq^ij?) is the running MS-mass, Cp = 4/3, ho = 25/3. 

5. For hght quarks the asymptotic forms of wave functions entering eq.(4- 
6) look as follows (see ch.9 and Appendix B in |9j). 

a) for the leading twist 2 wave functions: 

Pa{x) = Viix) = Vrix) = (pasyix) = 6X1X2 , (14) 

b) for the non-leading twist 3 wave functions: 

Pp{x) = l, V^ix) = ^[l + {x^-X2f], (15) 

Va{x) = Prix) = 6x1X2 . 

For heavy quarkonium the light-front IS'-Coulomb wave function can be 
taken as: ^ 

r , f kl + {l-^x,X2)M^' , 

0(x) r^jd^k^ n-, k^) - ^^^\ ^i_Cx7i-v-)] \ ■ ^'^^ 

Here: qs is the Bohr momentum and v = qb/Mq <^ 1 is the mean heavy 
quark velocity. 

^ It originates from the Schrodingcr equal-time wave function Schif) ^ exp(— gsr) 
^Sch(fc) ^ (I^P + "Zb)^^' supplemented with the proposed in jJU] and commonly used 
substitution ansatz: kj_ kj_, kz {xi — X2)Mq/2, Mq = {Mq"^ + k'\)/xiX2 ■ 
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For the leading twist charmonium wave functions we will use below a 
somewhat modified simple model form: ^ 



where v is now a parameter which has the meaning of the characteristic quark 
velocity in the bound state, while Co{v^) is the normalization constant. 

The wave function t;^) in eq.(17) interpolates in the simplest way 
between two extreme cases: very heavy quark with f — * in eq.(16), and 
very light quark with (formally) v ^ 1 in eq.(14). In the non-rclativistic case 
of very small v ^ 1 the wave function (j)o{x,v'^) is strongly peaked around 
the point Xi — X2 — 1/2, so that 0o(x,i;^ — > 0) — > S{xi — |). And clearly, 
a decreasing quark mass leads to larger v and to wider 0o(a;, f^). We take 
below = 0.30, as this value is commonly used in calculations for the IS* - 
charmonium. 

The wave functions 00 (a;, "i^^ = 0.30) fromeq.(17) and (pasyi^) = (po^x^v"^ = 
1) from eq.(14) are shown in Fig. 2 (the wave function (f)o{x,v'^ — 0.08) cor- 
responding to the scale of T(15') is also shown for comparison ^ ). 

Guided by the above examples, wc take for numerical calculations the 
following model wave functions of the 15'-charmonium: 



With = 0.30 this looks as follows (compare to eqs.(14,15), all wave 
functions are normalized: Jq dx(f)i{x, v'^) = 1) : 



It is inspired by the fact that after taking into account small relativistic corrections 
the Coulomb wave function of the 15 two-particle bound state behaves at small r typically 
as: 5'(r) ^ exp(— gsr). Here A = CqI"^ + O(w^). The constant cq is not universal 
and e.g. for the ^5*0 -states of hydrogen and positronium: Ahyd = -I- 0{v^), Apos = 
Av^ + 0{v'^). Therefore, this can be used only at really small v ^ 1, such that A ^ 1, 
and can't be taken literally for charmonium with ~ 0.3. So, wc have taken the simplest 
model form of eq.(17) having in mind that it behaves qualitatively in a right way and will 
be really apphed at ~ 0.3. 

^ It is seen from Fig.2 that even (j)o{x, = 0.08) is still far from the S - function. 




(17) 





(18) 
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Figure 2: The shape of the wave function 0o(a;,f^) for: = 1 - light quarks 
(asymptotic); = 0.30 - charmonium; = 0.08 - bottomonium. 
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a) for the leading twist 2 wave functions and Va{x): 

I \ 0.70 

Vt{x) = Vi^{x) = Pa{x) = Va{x) = 9.62 x^x^ . (19) 

\^ i — Z.O X\X2 ) 

b) for the non-leading twist 3 wave functions: 

/ X 0.70 

Pp(x) ^ 1.97 [ ^^^^ \ , 
^ ' 11-2.8x1X2/ 



/ X 0.70 

V.W = 1.67[l + fe-..n^-||l^ . (20) 



5. For numerical calculations we use the following parameters. As for 
Q;s(Ai^), we take the simplest form: (47r/feo) ln"^(/xV^^) with A = 200 MeV. 
It gives, in particular: as{iJ^ = M^) ~ 0.34 at the scale of the r-lepton 
mass, and this looks reasonable. For the MS'-mass of the c-quark we use 
the standard value: Mc = 1.2 GeV, and we take fp fv — 400 MeV on 
average ^ . 

With all this and the wave functions given above in eqs.(19,20), one ob- 
tains from eqs.(l,9,10): 

lo ~ 5.5 , Fvp{s = M^^^s)) ~ 3.5 • lO'^ GeV-\ 

a {e^e- ^ 7* ^ J/V'x + r}}j ~ 33 fb. (21) 

This agrees with the BELLE result (assuming that Br{r]c > 2chargcd) is 
not significantly less than unity). This fact can't be taken too literally, of 
course, due to dependence on the model form of used wave functions. Some 
check of the sensitivity of results to the wave function forms can be obtained 
by variation of the parameter in cq.(18). So, we repeated calculations with 
v"^ — (0.30 ± 0.05). Then the cross section changes only by ±7%. Therefore, 
the results are not very sensitive to reasonable variations of v^. 

^ The value of fv is known from the leptonie deeay of J/t/j: T{J/ip e^e^) = 
{IQtto^ /27){\Jv? /M^), \Jv\ ^ 410MeV. The one-loop eorreetions (- as/w) to the ratio 
T{r]c 2^)/r{J/ijj e+e~) indicate that fp is slightly smaller than fy. 
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Now, let us discuss some characteristic features of the above result for the 
integral /q in eq.(lO). The mean value of the renormalization factor Z^, see 
eqs.(12,13), is: (^m) — 0.74. This means that in Fig.l the typical value of the 
running c-quark mass in the quark propagator is: (Mc(/i^ = a^)) ~ 0.9 GeV, 
while (a^) ~ 33 GeV^ which is ~ 1.5 times smaller than the typical value 
~ ^ 46 GeV^ for the narrow wave functions. 

The value of (fc^) can be inferred either from the mean value of (Z^) ~ 
0.80, or from the coupling as{k'^) '■ {as{k'^)) — 0.263. Both give the same 
result and show that the mean virtuality of the gluon in Fig.l is: {k"^) ~ 
12GeVl This is ~ 2 times smaller than a typical rough estimate: (A;^) ~ 
g^/4 ~ 23 GeV^ for narrow S- like wave functions (and is not so far from the 
two-quark threshold (2M,)2 ^ 5 3 QeV^). ^° 

The smaller values of (fc^) and (cx^) is the main reason why the standard 
NRQCD-calculations underestimate the cross section considerably. In other 
words, the charm quark is not very heavy and its wave functions are not 
much like the (5-functions, although they are of course significantly narrower 
than similar wave functions of really light quarks (see Fig. 2). 

It is also of interest to make a comparison to the value of the cross section 
obtained from eqs. (1,9,10) in the limit which, in essence, corresponds to the 
approximations of NRQCD, see ^]). For this, one has to replace in eq.(lO): 
a) all wave functions, except for Va{x) , hj Vi = Pi = 6{x — 1/2); b) to omit 
the term with Va(x); b) to replace all by 1. In this case: 

~ [1 + (2MQ/go)2 | ~ 26 GeV^ a^i^) ~ 0.23 . (22) 

Substituting all this into eqs. (1,9, 10), one obtains: 

/o~1.6, Fv^p ~ 1.1 ■ lO'^GeV^S (T~3fb, (23) 
which is essentially smaller than in eq.(21). 

On the whole, as was argued above, the difficulties with explaining the 
BELLE result for (j(e"'"e' J/ip + rjc) are not really the difficulties of QCD, 

^°It is worth noting also that while /°(/i*) in eq.(7) is very small and negative (e.g., 
using for the c-quark pole mass the value M* = 1.65 GeV obtained in 11 from the D- 
meson semileptonic width, one obtains: /"(/i*) — —0.08 ), the coupling /°(fc^) from eq.(8) 
becomes positive and much larger, (/^(fc^)) — +0.26, so that the contribution of the term 
with the wave function Va{x) in eq.(lO) is ~ 15% of the whole integral Iq- 

In the opposite limit when all IS* - charmonium wave functions are taken as the 
asymptotic ones, see eqs. (14, 15), the cross section will he a ~ 70 fb. 
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but are rather due to a poor approximation of the real dynamics of c-quarks 
by NRQCD. Within the approach described in this paper (which we consider 
as more reahstic), the BELLE results look rather natural. 

We hope that subsequent experimental and theoretical efforts in this field 
will help to elucidate dynamics and properties of various charmonium states. 
In particular, the use of the wave functions Vi{x) given in eqs. (19,20) instead 
oi 6{x — 0.5) from NRQCD, will enhance the calculated cross sections of 
inclusive direct production of J/ip. 

We are grateful to A.I. Milstein for explaining us the properties of the 
positronium wave functions and a useful discussion. 

The work of V.L. Chernyak was supported in part by the RFFI grant 
03-02-16348-a. 
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